Grain boundaries in vortex matter 
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We explore the statistical properties of grain boundaries in the vortex poly crystalline phase of type 
II superconductors. Treating grain boundaries as arrays of dislocations interacting through linear 
elasticity, we show that self-interaction of a deformed grain boundary is equivalent to a non-local 
long-range surface tension. This affects the pinning properties of grain boundaries, that are found 
to be less rough than isolated dislocations. The presence of grain boundaries has an important effect 
on the transport properties of type II superconductors as we show by numerical simulations: our 
results indicate that the critical current is higher for a vortex polycrystal than for a regular vortex 
lattice. Finally, we discuss the possible role of grain boundaries in vortex lattice melting. Through a 
phenomenological theory we show that melting can be preceded by an intermediate polycrystalline 
phase. 



I. INTRODUCTION 

The phase diagram of high temperature superconduc- 
tors is an object of active investigation in condensed mat- 
ter physics. Most high Tc materials behave in a magnetic 
field as type II superconductors, with further complica- 
tions due to the broader phase space — in terms of tem- 
perature T and field H — in comparison to conventional 
type II superconductorsiSi^. This leads to several pos- 
sibilities for the mixed phase, where magnetic flux pen- 
etration is incomplete. As first discussed by Abrikosov 
for conventional superconductora^i, flux is quantized and 
carried by vortex lines which are arranged in the form 
of a lattice. As in conventional matter strong enough 
fluctuations destroy long range order: when temperature 
is raised the vortex lattice melts into a vortex liquid^'S'i. 
Fluctuations are also provided by quenched disorder that 
is intrinsically present in these materials, leading to com- 
plex glassy phasesSiSiifliiiiiSiii. 

While several experimental methods have been used 
to investigate vortex matter, a direct image of the geo- 
metrical and topological properties of the vortices can be 
obtained by the Bitter decoration technique^'*. Its ap- 
plication to conventional superconductors provided the 
flrst direct proof of the vortex latticei^ predicted by 
Abrikosovi. The observed lattice contains, however, 
topological defects, such as dislocations and grain bound- 
aries. These last extended defects are the signature 
of a vortex polycrystal with crystalline grains of dif- 
ferent orientationsi^ii^. Vortex polycrystals have been 
observed, after field cooling, in various superconducting 
materials such as NbMcl''■l^ NbSei"-^*-^^-^", BSCCO^i 
and The grain size is typically found to grow 

with applied magnetic fieldiSiii. Moreover, two-sided 
decoration experiments show that the grain boundaries 
thread the sample from top to bottonii£*i^, i.e., one ob- 
serves a columnar grain structure. Despite the wealth 



of experimental observations, there is no detailed theory 
accounting for the formation of vortex polycrystals. 

The behavior of vortex matter in presence of disor- 
der represents a formidable theoretical problem that has 
still not been completely solved. While early theoret- 
ical considerations seemed to imply that even a small 
amount of disorder would lead to the loss of long-range 
order— and to the formation of an amorphous vortex 
glass phas©ii, it is now accepted that at low disorder 
vortices arrange into a topologically ordered phase: the 
Bragg glassi^iSi. The existence of this phase, character- 
ized by logarithmically growing correlations, slow relax- 
ation, and other glassy features, has been now experimen- 
tally confirmed^^. At high enough disorder, the Bragg 
glass phase is found to be unstable against dislocation 
proliferation and one may expect the transition into an 
amorphous vortex glasa^SiSLSS. The precise nature of this 
transition and, more generally, the mechanism underly- 
ing vortex lattice melting is still under debate. Typical 
melting theories are based on variants of the Lindemann 
criterion with disordes^S, or involve dislocation prolifer- 
ation mechanisms^Si. 

The properties of dislocations in the vortex lat- 
tice have been the object of extensive theoretical 
investigations^iiSSiSi^, but grain boundaries are less 
studied although they are often observed in numerical 
simulations^Si^Si^l. For instance, the vortex plastic flow 
in the Corbino disk geometry is characterized by radial 
grain boundaries sliding in the tangential directioni^ . In 
addition, recent numerical simulations indicate the pres- 
ence of an intermediate polycrystalline phase before the 
melting transitionSS*^. This behavior was observed us- 
ing different numerical methods in two dimensions^ and 
in presence of columnar disordep3i. This suggests that, 
in some conditions, grain boundaries may play a role in 
the melting process, as in the theory of grain boundary 
induced melting of two dimensional crystals^. 
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Here we analyze the properties of grain boundaries in 
vortex matter describing the fluctuations induced by dis- 
order, stress or temperature. A grain boundary can be 
considered as an array of dislocations, whose dynamics is 
ruled by internal stresses. While ideally a grain bound- 
ary minimizes its energy by remaining flat, the action of 
external perturbations leads to deformations that can be 
described by the theory of elasticity^^. We compute the 
self-interaction of a deformed grain boundary extending 
the results obtained for isotropic elasticitj*^ to the case 
of the vortex lattice. Grain boundaries are much stiffer 
than isolated dislocations, possessing a non-local long- 
range surface tension and, in presence of disorder, they 
are expected to be less rough than isolated dislocations. 
We estimate the grain boundary roughness exponent us- 
ing the random stress model introduced in Ref. 34 for 
vortex dislocations. Using scaling arguments, we also 
derive the creep law for thermal activated motion and 
discuss disorder arrested grain growths. 

The critical current is an important property of type 
II superconductors, since it represents the current be- 
low which vortices are pinned and the material conducts 
without resistance. It is thus interesting to understand 
how the topological properties of vortex matter influence 
its behavior. We use numerical simulations of interact- 
ing vortices to quantify the effect of grain boundaries 
on the critical current. We obtain a polycrystalline vor- 
tex structure by relaxing at zero temperature a random 
initial vortex arrangement. This process simulates a typ- 
ical field cooling experiment in which the temperature is 
rapidly decreased from above Tc in presence of a field. 
The system moves rapidly towards lower energy config- 
urations corresponding to zero temperature and thermal 
effects can thus be disregarded. In this case magnetic 
flux is present in the material as it enters the supercon- 
ducting phase and vortices are initially disordered. Once 
grain growth has stopped, we simulate the effect of an 
external current flowing through the sample by apply- 
ing a constant Lorenz force. The critical current is then 
defined as the current at which vortices start to move 
steadily. By repeating the simulations for different val- 
ues of the vortex number, representing the effect of var- 
ious magnetic field intensities, we show that the critical 
current for a polycrystal is always larger than the one ob- 
tained for a perfect lattice. This reflects the fact that a 
polycrystalline assembly is more effectively pinned than 
a perfect lattice because it can accommodate better in 
the disordered landscape. In addition, we find that the 
corresponding IV curve is hysteretic upon ramping up 
and down the current. This result can explain the differ- 
ence in transport properties between field cooled and zero 
field cooled samples and the related hysteresis commonly 
measured experimentall}ii2i^i*^2iiM4iiS4&. 

Finally, we discuss the possible role of grain bound- 
ary in vortex lattice melting by constructing a free en- 
ergy functional for the grain boundary density along the 
lines of Ref. 's^. We derive the contribution due to grain 
boundary fluctuations and junction formation and show 



the presence of a polycrystalline phase with a finite grain 
boundary density. As temperature is increased the grain 
boundary density increases and the system melts. The 
theory predicts a value for the melting temperature that 
is quite similar to that obtained in Ref. ,30 considering 
only the contribution of isolated dislocations. We no- 
tice that the presence of an intermediate polycrystalline 
phase in vortex lattice melting was recently proposed on 
a phenomenological basis in Ref.**?! and experimental ev- 
idence was reported for Lai.gSrg.iCuO^^. 

The paper is organized as follows: in section^]we com- 
pute the self-energy of a deformed vortex grain boundary. 
In section IIIII we analyze the interaction between grain 
boundaries and disorder, as well as its relevance for pin- 
ning, creep, and grain growth. Section IIVI reports the 
results of numerical simulations of interacting vortices 
where we discuss the effect of grain boundaries on the 
critical current. In section we discuss the role of grain 
boundaries in the melting process through a phenomeno- 
logical theory. Section is devoted to conclusions. Fi- 
nally, the appendix reports details of the derivation re- 
ported in section [V] 



II. ELASTICITY OF GRAIN BOUNDARIES 

A simplified but rather effective description of the vor- 
tex lattice is provided by its representation as an elastic 
crystal of flux lines. At large enough distances, the elas- 
tic energy of the vortex lattice can be expressed in terms 
of the vortex displacement field u as follows 

n^^Jd'^r [C66(VU)2 + (cii - C66)(V • U)2 +C44(a,u)2] , 

(1) 

where cn, C44, cge are the local elastic moduli, and 
the magnetic induction B is parallel to the z direction. 
Within this representation, we shall introduce an ideal 
low angle grain boundary as an infinite periodic array of 
straight dislocations in the vortex lattice oriented along 
the z axis, spatially arranged along the y axis with an ar- 
ray spacing equal to D, and with Burgers vectors b point- 
ing along the x direction (i.e. edge dislocations). The 
wandering of the i-th dislocation line can be schematized 
through the vector Ri(z) = {Xi + Xi{z),iD), assum- 
ing that all displacements take place within glide planes, 
i.e. the xz plane, so that Xi + Xi{z) plays the role of 
the displacement field of the grain boundary as well. Xi 
is a constant term and deals with rigid displacements 
of the dislocation lines. Its contribution to the elastic 
Hamiltonian is known since it is the same as for straight 
dislocations in isotropic lattices'^^. In the following this 
contribution will be referred to as Tip- 

Defining = the vector u can be decomposed 

as u(r) = u''(r_L, u?(r_L— Ri(z), z), where u?(r_L — 

Ri(z),z) is the singular solution of the two-dimensional 
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problem for each value of z 

r ceeV^uf + (cn - C66)V(V -0=0 . 

while u'"(r^, z) is the regular part of the solution due to 
the interplane couplings along z. 

Minimizing Eq. with respect to u and imposing the 
first expression of Eqs. (jSJ we find the differential equa- 
tion 

C66V2u'' + (cil-C66)V(V • U'')+C44<9fu'' = -C44<9^ ^ U,^ 

i 

(3) 

where the field u| on the right-hand side term of the 
equation is known from elasticity theory as the displace- 
ment field generated by a point edge dislocation at Ri(2;). 
Performing a first order expansion in the displacement 



Xi^Xi{z)^ the derivative removes any dependence on the 
constant part X^, and we can rewrite Eq.® in Fourier 
space as follows 

^2 

cgeg^u'^ + (cii - C66)q(q ■ u'^) + caklu'' = c^^^A, (4) 

q 

where q = {kx,ky) and 

A _^\^^^kyY^x (k ] ( (1 -^)cos20] \ , . 

A-2^e ^"(^-M A:,[r + (f-r)cos20] j 

with r = cge/cii, cos0 = k^jk and sin(/) — ky/k. 

A can be decomposed in its longitudinal and transverse 
components = q(q • A)/q^ and At = A — Al. The 
Hamiltonian ^ thus becomes 



(6) 



where we have neglected constant terms and defined 



M(g,</), k,) 



C66 cos-^ 20 Ciir^ sin^ 20 



(7) 



Defining Xn{kz) = / e~''3a"^ X{Qy, k^), where the integral is restricted to the first Brillouin zone (BZ), we 

Jbz 27r 

get 



^^r^'^T^J^^^J ^my + Gy,k,)XiQy,k,)Xi^Qy 



(8) 



where we have introduced the interaction kernel 



/ + 00 
M{k:„Qy+Gy,k,) 



with 



M[kx, ky, kz) 



4r^ 



{kl + kinki + ki + {kl + kiYiki + ki + ^^^ki) ■ 

Solving the integral in Eq.© leads to 
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= n" + 


7r&2 


2D^ C44 
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"(2A^4 


Cee ^ J 


k1 




+ ^k1 

Cee ^ 



Cll 



C66 Cll 



with ky — Qy + Gy, S> f and ^ 1. 

C66 .Cll 

Moreover, keeping the leading term of the righthand side in Ea. (|f l|l we get 



(9) 



(10) 



(11) 
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(12) 
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It is a common procedure to rescale the y coordinate 
by a factor — 1 / — in order to get an isotropic reference 

2 V C66 

frame. The elastic Hamihonian thus becomes 



H = K 



TTb^ 

2L>2 



E 



(2^ 



|k|X(k)X(-k), 



(13) 



being k = {ky, k^) and K = y'ciJcee- 

In this Hmit, the same result predicted by the isotropic 
theorji^S is thus obtained. The nonlocal character of the 
elastic kernel (oc k) manifests that long range interactions 
between dislocations stiffen the grain boundary, and that 
a surface tension approximation is not suitable for a cor- 
rect description of its elastic properties. 



III. INTERACTION BETWEEN GRAIN 
BOUNDARIES AND DISORDER 

A. Random stresses 

Point defects such as vacancies or interstitials in the 
underlying crystalline structure of the superconducting 
material, and/or substitutional impurities, etc., act as 
pinning centers for the magnetic vortices. For weak pin- 
ning forces, disorder can be theoretically described by 
a random pinning potential acting directly on flux lines. 
The distortions generated in the vortex lattice as well 
as the occurrence of depinning under an applied current 
have been intensively studied over the last decades (for 
a review see Ref. Q). 

Here, we are instead concerned with the behavior of 
grain boundaries in presence of disorder. The disorder in- 
duced vortex lattice displacement field gives rise to shear 
elastic stresses, which, in turn, generate Peach-Koehler 
forces on the vortex lattice dislocations^. In other words, 
as the final consequence of these disorder-induced dis- 
tortions, there is an effective pinning stress field crij(r) 
acting as well on vortex dislocations (and therefore on 
grain boundaries). The statistical properties of the ran- 
dom stress field has been analyzed in Ref. in the case 
of vortex dislocations. In the following, we recall their 
derivation and adapt it to the case of grain boundaries. 

On short length scales, where vortex displacements 
u(r) are smaller than the coherence length f (the so 
called Larkin regime^), a perturbative calculation can 
be performed. As discussed in Ref. Isol for grain bound- 
aries it is necessary to consider larger scales, ^ < u < a, 
where vortices are well described by a Random Manifold 



(RM) modeliJ^ in which flux lines are subject to an un- 
correlated pinning potential. In this case, the relative 
displacements correlation function is 



i?,,(r - r') = [u,ir) - ^.(r')]K(r) - ",(r')] 



(14) 



Here Ra is the crossover length, also known as positional 
correlation length, at which average vortex displacements 
are of the order of a. The roughness exponent can be 
estimated as (^rm ~ 1/5. 

On scales larger than Ra, vortex displacements are of 
the order of a and the periodicity of the lattice comes 
into play^'^. Displacements are shown to grow loga- 
rithmically, with correlations of the form B{r — r') ~ 



-I'ln^"'^ 

TT 



Ra 



and topological defects are absent. 



This quasi-ordered phase is knwon as the Bragg glass 
(BrG)i3. 

The defect-free regions discussed above act on vor- 
tex lattice dislocations through a Peach-Koehler stress 
field^. Statistical properties of this stress field can be 
obtained from the correlator Bij{r — r') applying lin- 
ear elasticity theory. In particular, the stress correlator 
Sxy{r - r') = axy{r)axy{r') wiU read 



S,y{r - r') = (i^V2) [dA'Byyir - r')+ 
dydy'Bxx{r - r') + 2dxdy'Byx{r - r')] . 



(15) 



Replacing previous expressions of Bij (r — r') we easily 
obtain the stress fluctuations over a distance R 



i?2 



{R/Raf^'"'' R<Ra 
1 R> Ra 



(16) 



where the first case applies to the RM description, while 
the second corresponds to the BrG regime. The effect of 
this random stress on isolated dislocations was studied 
in Ref. 34 where several differences with respect to the 
case of vortex lines were pointed out. Here we consider 
the behavior of grain boundaries, expecting substantial 
novel features arising from long range interactions be- 
tween grain boundary dislocations. 

The Hamiltonian of a grain boundary in presence of 
disorder can be written as 



Hd —'H + Hpin, 



(17) 



with TL being the elastic term calculated above and TLpm 
the pinning term given by 



'Hpin = ^ j dzXi{z)h<7xy[Xi{z),iD, 



(18) 
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Although there is no explicit expression for 7i„m, it is possible to derive its fluctuations over a distance L as 



J 



Ep^n^b'^y^ dz' X,{z)X^, {z') a^y[X,{z), iD, z]a^y[X^, {z'),i'D, z']. 

i ,/ "'0 Jo 



(19) 



r 



Taking the continuum limit of the sum and integrating 
for both the RM and the BrG regimes, the typical pin- 
ning energy when displacing a grain boundary segment of 

length L by an amount Xf ^ uqe will be thus given 

by 



Kab\ 9 9 



In- 



(RM) 
(BrG) 



(20) 



A dimensional estimate of the elastic cost of fluctuations 
of a grain boundary fraction of linear dimension L has 
the form 



E, 



(21) 



Since we are dealing with static properties of the sys- 
tem, we can impose equilibrium conditions balancing E 
and Epin, that is, equating the elastic cost of fluctuations 
and the energy gain due to the interaction with disorder. 
Defining the roughness exponent of a grain boundary C,qb 



from u^g 



•,GB 



i (RM) 
log^/^ (BrG) 



(22) 



The long-range stiffness of a grain boundary reduces the 
values of roughness exponents in comparison with the 
case of isolated dislocations^. 



B. Depinning and creep 

So far we have not considered the effect of driving 
forces on the dislocation arrangement. Driving forces for 
grain boundary motion can be externally induced by a 
current flowing in the superconductor or internally gen- 
erated by the ordering process during grain growth'"^ ^. In 
both cases, the presence of a driving shear stress a gives 
rise to a Peach-Koehler force per unit length of the form 
Fdrive = crb acting on each dislocation along the grain 
boundary fraction considered or, in other words, to a to- 
tal driving force per unit length equal to Fdrive = abL/D. 

At low stress grain boundaries are pinned. One can 
estimate the depinning stress from conventional scaling 
arguments. The energy associated to the driving force 



acting on a low-angle grain boundary segment of length 
L and displaced by an amount uqb is given by 



E, 



drivi 



■XL) 



E 



dzF'^„^^{z)uGB{yz,z) 



abL^ 



D 



-UGB- 



(23) 

The depinning stress can be obtained comparing this 
driving term with the pinning energy reported in Eq. 1201 
The relevant scale to consider is due to the interplay be- 
tween elasticity and disorder and results from the min- 
imization of Eel + Epin for displacements of the order 
of UGB ~ a ~ 6, corresponding to the dislocation core. 
A similar approach is followed in the case of vortices^, 
which are pinned for displacements of the order of ^, the 
size of the vortex core and, hence, the relevant scale for 
the interaction with impurities. In our case, we obtain 
the Larkin length as Lp ~ {b/D)^Ra, which is typically 
smaller than Ra- The depinning stress is then identified 
as the stress necessary to depin a section of dimension 



C7e ~ Kby{DLp) = KD^ib'^Ra). 



(24) 



For low values of the stress (cr ^ CTc), the response of 
a grain boundary is mainly due to thermally activated 
motion in a disordered environment^. In this case, we 
expect a highly non-linear creep motion with an aver- 
age velocity v ~ exp[~C{(Jc/(j)^ /T], where C is a con- 
stant, and ^ is the creep exponent that quantifies the 
divergence of the energy barriers U{a) ~ separat- 
ing metastable states. An estimation of the exponent 
fi for a grain boundary can be obtained from a simple 
dimensional scaling argument, which is confirmed by a 
more rigorous renormalization group analysis. The typi- 
cal energy barrier for a grain boundary section of length 
L is of the order of U (L) ^ E^+'^C.gb ^ -^Jigre we have 
used ug b ~ L^*^^ . In presence of an applied stress 
cr, we can compute the typical grain boundary length 
L{a) involved in thermal activated motion minimizing 
U{L) + Ed„ve{L). The result yields L{a) - a^/iCcB-i)^ 
Using this length, we obtain that the typical energy bar- 
rier depends on the stress as 

^(1+2Cgb)/(Cgb-1) 



(25) 



implying that fi = 2Cgb/(2-Cgb)- For the RM and BrG 
regimes the exponents are given by 



I (RM) 
1 (BrG) 



(26) 
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exponent 


length scale 


isolated dislocation 


2D bundle 


grain boundary 


c 


RM 


15/13 


5/13 


1/5 


c 


BrG 


1 - log2/3 


1/3 


log^/^ 




RM 


17/11 


10/21 


7/4 




BrG 


1 


2/5 


1 



TABLE I: Comparison between roughness and creep exponents calculated for isolated dislocations, 2D dislocation bundles—, 
and low angle grain boundaries, taking into account non-local effects proven in Section HTl 



Now these exponents are larger than their counterparts 
calculated for isolated dislocations. In other words, the 
formation of grain boundaries affects vortex dynamics 
lowering ordinary creep rates. On table llTr^ all previous 
results are summarized and compared to estimates for 
different dislocation arrays. 



An explicit expression for the grain size can be obtained 
inserting ac^ reported in Eq. 1241 for the weak pinning 
regime, yielding Rg ex Ra- A similar calculation can be 
performed in the strong pinning regime^ and the results 
appear to be in good agreement with experiments mea- 
suring average grain sizes in NbMoi&. 



C. Grain growth 



IV. THE EFFECT OF GRAIN BOUNDARIES 
ON THE CRITICAL CURRENT 



In a field cooling experiment, magnetic flux is already 
present in the sample as it is quenched in the mixed 
superconducting phase. It is thus reasonable to expect 
that vortices are originally disordered and that, due to 
their mutual interactions, undergo a local ordering pro- 
cess. Along this process, many dislocations annihilate, 
and most of the remaining dislocations arrange them- 
selves into grain boundaries with various orientations. 
The growth of crystalline vortex grains is due to the 
motion of these separating boundaries. The resulting 
polycrystalline structure has been indeed observed exper- 
imentally by means of Bitter decorations of both highSLS^ 
and low j'jiiigiiS superconducting samples. The effect of 
quenched disorder is to pin the grain boundaries, hinder- 
ing the growth process. Thus to understand the proper- 
ties of vortex polycrystals, it is important to analyze the 
dynamics of grain boundaries in vortex matter as they 
interact with disorder—. 

Grain growth is driven by a reduction in energy: For an 
average grain size R and straight grain boundaries, the 
characteristic energy stored per unit volume in the form 
of grain boundary dislocations is of the order of Fo/i?, 
where Fq is the energy per unit area of a grain bound- 
ary. Hence, the energy gain achieved by increasing the 
grain size by dR is T^/R^dR. Physically, the removal of 
grain boundary dislocations occurs through the motion 
of junction points in the grain boundary network. As 
junction points must drag the connecting boundary with 
them, which may be pinned by disorder, motion can only 
occur if the energy gain at least matches the dissipative 
work which has to be done against the pinning forces. 
The dissipative work per unit volume expended in mov- 
ing all grain boundaries by dR is aJ)/{DR)dR, where <Tc 
is the pinning force per unit area. Balancing against the 
energy gain yields the limit grain size 



Rn 



DTo 
bar 



(27) 



Numerical simulations of interacting vortices in two di- 
mensions (2D) allow to verify and to keep track of the 
ordering process and the grain formation after a rapid 
field cooling of the vortex system in thin superconduct- 
ing films. To this end, we consider a square 2D supercon- 
ducting cross-section of linear dimension L perpendicu- 
lar to the external magnetic field B along the z direction, 
where we locate a set of Ny rigid vortices (for most of the 
results presented here, we have considered values of N^j 
ranging from 516 to 4128). The dynamics of each vortex 
line i at position can be described by an overdamped 
equation of motion of the form 

Tdv,/dt = ^f™(r,-r,)+^f,p(r,-rp+fi(r,), (28) 



where F is an effective viscosity for vortex flow. The 
first term on the right hand side of this equation fol- 
lows from the fact that a pair of vortices interact with 
each other via a long-range force fui,(r) — AKi{\y\/ X)f ^ 
where A = $Q/(87r^A'^), $0 is the quantized flux car- 
ried by the vortices, A is the London penetration length, 
and Ki is a first order modified Bessel function^. Dis- 
tances are always measured in units of A. The second 
contribution reflects the attractive interaction forces be- 
tween vortex lines and quenched in point defects such 
as oxygen vacancies or other impurities in the material. 
These pinning centers are randomly located at positions 
(i = 1, . . . , Np) within the simulation box, and exert 
pinning forces according to a Gaussian potential of the 
form V{t — v'P) = Vbexp[— (r — r^)^/^^], whose ampli- 
tude and standard deviation are Vq and ^, the character- 
istic coherence length of the superconductor, respectively 
(The usual number of pinning centers Np considered is 
4128, and we have chosen ^ = 0.2A, characteristic of low 
temperature superconductors such as NbSe or NbMo.). 
Finally, if an external current J(r) is eventually applied 



FIG. 1: Relaxation of the topological defect structure from a 
simulation of Nv — 4128 interacting vortices after a sudden 
field cooling from a disordered vortex state in a simulation 
cell of linear size L = 36A. The colored five/seven-fold co- 
ordinated vortices (filled circles) indicate dislocations in the 
vortex lattice. The final configuration (snapshot (d)) is com- 
pletely pinned by disorder. There one can observe a polycrys- 
talline structure with most dislocations arranged into grain 
boundaries. 



FIG. 2: Pinned vortex structure for different values of the 
magnetic field: (a) = 1460, (b) A''^ = 2064, (c) iV„ = 
2919, (d) = 4128, after a sudden field cooling from a 
disordered vortex state in a simulation cell of linear size L = 
36A. The colored five/seven-fold coordinated vortices (filled 
circles) indicate dislocations in the vortex lattice. The average 
grain size in the resulting polycrystalline structure seems to 
grow with the intensity of the average magnetic field inside 
the cell. 



to the sample, it generates a Lorentz-like force acting on 
the vortices fL(r) = <I>oJ(r) x z/c, where c is the speed of 
Hght. These coupled equations of motion 1^5)1 are numer- 
ically solved with an adaptive step-size fifth order Runge- 
Kutta algorithm, imposing periodic boundary conditions 
in both directions. 

We first consider the relaxation dynamics of the vortex 
hues in the absence of driving currents. Moreover, in the 
present analysis we completely disregard thermal effects, 
that is, we mimic the dynamics of the vortex system af- 
ter a sudden quench of the superconducting sample from 
high temperatures (or equivalently, random vortex con- 
figurations) towards the lower energy states correspond- 
ing to zero temperature. After a transient regime, the 
dynamics stops due to disorder. We analyze the result- 
ing spatial configuration of flux lines by means of Delau- 
nay triangulations. A pair of a five-fold and a seven-fold 
neighboring vortices correspond to a dislocation in the 
vortex lattice. In the course of the simulations, the num- 
ber of five/seven-fold coordinated vortices is the same, in- 
dicating that during the relaxation process no other topo- 
logical defects such as disclinations appear to be present 
in the lattice. 

In Fig. n] we report a series of snapshots illustrating 
that the gradual ordering process involves the arrange- 
ment of dislocations in grain boundaries. The formation 
of these walls of dislocations screens out the long range 
elastic stress and strain fields otherwise created by dis- 
locations in the lattice and, at the same time, they ren- 
der a polycrystalline structure of the vortex array. This 



polycrystalline structure evolves in time until the residual 
stresses accumulated in the distorted vortex lattice drop 
down below the critical value CTc. At this point, grain 
boundaries get pinned by disorder limiting the average 
grain size (see Fig.^d)). Moreover, the limit grain size 
Rg/a appears to increase with magnetic field B cx Ny 
(see Fig. EJ, in qualitative agreement with experimen- 
tal results and the theoretical predictions reported in 
Ref.E 

Different experimental, or simulation, protocols will 
certainly influence the relaxation dynamics and the re- 
sulting metastable configurations of trapped disloca- 
tions and grain boundaries. Metastability and history- 
dependent features have been long recognized in driven 
vortex lattices^. We have considered a field-cooling pro- 
cedure since most of the Bitter decoration experiments 
are performed in a similar manneriiSiSSiS, and can thus 
be well described by the current simulations. Neverthe- 
less, other numerical protocols can also be devised, as for 
instance the one recently proposed in Ref. |3^ to exam- 
ine the vortex topology across the so-called peak-effect, 
that are better suited to reproduce diverse experimental 
conditions. 

Next, we study the behavior of the critical current 
Jc{B) for these 2D vortex polycrystals by means of nu- 
merical simulations. An externally applied current may 
induce the annealing of the metastable configurations (at 
least, to a certain degree that obviously depends on its 
intensity) present in Figs. ^ and [3 This is indeed ob- 
served in our numerical simulations, where we can as 
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FIG. 3: Vortex trajectories between two pinned configura- 
tions obtained after the application of a small driving current 
below the threshold value Jc{B). Small and heterogeneously 
distributed displacements of the vortex positions are observed 
in both the parallel and perpendicular direction to the applied 
force /£. The number of vortices in the simulation cell of lin- 
ear size L = 36A is Nv = 2919. The number of pinning points 
Np = 4128. 



well identify the critical current Jc below which the av- 
erage motion of the vortex lattice eventually ceases after 
a rich initial transient of plastic flow. As one can ob- 
serve, for instance, in Fig. 13 a small current below the 
threshold value Jc{B) gives rise to non-trivial (i.e. not 
just a slight drift along the force direction /£) changes of 
the displacement field Un of the vortex lattice. The vor- 
tex displacements are heterogeneously distributed and a 
small component perpendicular to the force direction can 
be observed. This in turn, implies changes of the elas- 
tic shear stress distribution responsible for the Peach- 
Koehler forces acting on grain boundary dislocations 
that, as a consequence, may move and rearrange in re- 
sponse to the new force field. 

We have determined the dependence of the critical cur- 
rent on the magnetic field by carrying out simulations 
for different densities of vortices in the simulation cell. 
Moreover, we have compared these results with those 
corresponding to a completely different initial state: a 
perfect single crystal configuration with similar densities. 
Our results are summarized in Fig. ^ The qualitative 
and quantitative differences between the two curves rep- 
resented in the figure are due to the presence of grain 
boundaries. The presence of these topological defects 
in the vortex configuration enhances the critical current 
needed to give rise to a steady regime of plastic flux flow, 
that in this case, appears to be controlled by grain bound- 
ary motion. The plastic deformation of crystals is usually 
mediated by the nucleation and motion of dislocations^. 
Another possible mechanism for plastic flow is the glide 
motion of grain boundaries which, as in this case, can 
be the most relevant mechanism when the grain sizes are 




1000 2000 3000 4000 5000 

FIG. 4: The critical current Jc as a function of the number 
of vortices A''„ in the simulation cell. The number of pinning 
points Np — 4128, the cell size L — 36A, and the Ginzburg- 
Landau parameter considered is k = 5. The upper line shows 
the results obtained starting from initial field-cooled configu- 
rations containing grain boundaries (GB), whereas the lower 
curve shows the numerical results obtained from perfect crys- 
talline initial configurations. 



limited and there is a high fraction of grain boundary 
atoms. According to our numerical results, grain bound- 
aries are more efficiently pinned by disorder, in agreement 
with the general expectation that grains adjust better to 
the disordered landscape than a perfectly ordered lattice. 
In both cases, we observe the decrease of Jc with an in- 
creasing density of vortices until this reaches a plateau 
for the largest number of vortices considered. 

It is also worth noting that we have not considered the 
renormalization of either the penetration length A or the 
coherence length ^ of the superconductor with the inten- 
sity of the magnetic field B. Within a mean-field sce- 
nario, these parameters should diverge as the magnetic 
field approaches the upper critical field = ^o/^ir^^. 
An estimation of the reduced field values we are dealing 
with in the simulations yields B/Bc^ = 2TTNy/ {k^V^) ^ 
(0.1 — 0.8). This means that the renormalization of A 
and ^ will be especially relevant for the last point of the 
curves in Fig. 01 Recent simulations of similar vortex 
lattices in 2D'^^ show that indeed such a field renormal- 
ization could be responsible for a sudden increase of the 
critical current close to the upper critical field Bc2 ■ 

On an experimental ground, our results match, at least 
on a qualitative basis, the behavior exhibited by vortex 
matter in critical current measurements at low magnetic 
fields. As stated above, grain boundaries are commonly 
observed in field-cooled (FC) samples. On the other 
hand, ordered vortex crystals can be obtained in zero field 
cooling (ZFC) experiments, i.e. applying a magnetic field 
only after temperature has been lowered to the expected 
value^"i^^i'^^i'^3,44^45,46^ g^^^g usually charac- 
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FIG. 5: The steady state average velocity of the vortices as 
a function of the apphed current J. The current is ramped 
up (and down) in steps and is kept constant after each step 
until the system reaches a steady state. The arrows indicate 
the direction of the ramp. The number of vortices is = 
2064, the number of pinning points A'p = 4128, the cell size 
L = 36A, and the Ginzburg-Landau parameter considered is 
K = 5. 



terized by a higher critical current and has been proven 
to be metastabl9^2ii4. These aspects result in a peculiar 
hysteretic behavior commonly observed in critical current 
measurementa^Siii and I — V characteristics^i*^. In our 
numerical analysis, the evaluation of critical currents in 
perfect vortex crystals (lower line in Fig.^l fairly mimics 
phenomenology of ZFC measurements, while results for 
the grain boundary model (upper line in Fig. ^ can be 
interpreted as a simulation of FC response. Hysteresis 
is in fact reproduced by our simulations when we start 
from the polycrystalline state. As shown in Fig. |S1 when 
the current J is ramped up vortices start to move at a 
current Jd , with a velocity that then increases with the 
current. If the current is ramped down from the moving 
state, vortices get pinned at a lower value of the cur- 
rent Jc2 corresponding to the critical current measured 
for a perfect crystal upon ramping up the field. Notice 
the similarity with the experimental results of Refs4i*^. 
Once more, we should underline how these results hold 
only for low values of the applied field. As the magnetic 
induction approaches its critical value, a sudden increase 
in measured critical currents is observed in both the ZFC 
and the FC experimental setwp—^. 



determines the transition to a liquid phase^'^'"^, while the 
effects of disorder associated with high magnetic fields 
are responsible for the insurgence of a glassy phas o^i^i^" . 
A deep theoretical understanding of such transition phe- 
nomena, accounting for their microscopic origin, has not 
been achieved yet. Nonetheless it has been shown that for 
strong enough disorder the Bragg glass phase is unstable 
against dislocation formation^^'^^''^^. This suggests that 
the melting process could be ruled by topological defects 
(as discussed in Ref. IboI) in analogy with two dimensional 
theories of crystal melting. Here we discuss the possibil- 
ity that, under the effects of fluctuations, dislocations 
unbind and rearrange in grain boundaries giving rise to 
a polycrystalline structure^. In this framework, the vor- 
tex polycrystal can be seen as an intermediate stage in a 
process that ends in the amorphous or liquid phases. 

Our purpose is to study the quasiequilibrium proper- 
ties of such polycrystalline stage, using the elastic proper- 
ties of grain boundaries in a vortex lattice derived in Sec- 
tion II. The main goal of our analysis is to write the free 
energy density / of the system as a function of different 
lattice arrangements in configuration space. A minimum 
in free energy for a polycrystalline configuration in prox- 
imity of the melting line would corroborate the hypoth- 
esis of a grain boundary mediated transition. For our 
purposes, we parametrize configuration space in terms of 
linear grain boundary density n, meaning that a n ^ 
configuration corresponds to an ordered (grain bound- 
ary free) vortex lattice. Our consideration focus on the 
thermally induced melting transition and the effects of 
impurities are neglected. 

We consider arrays of edge dislocations, parallel to the 
z axis and arranged in low angle grain boundaries. As in 
the case of grain growth, all Burgers vectors are in the xy 
plane, corresponding to a columnar grain structure. Fol- 
lowing the aforementioned ideas'^*, we can introduce the 
linear concentration of grain boundaries n and in the low 
density limit we can expand the free energy functional 
(per unit volume) in powers of n as 



/(n) = (70 + 7t) 71 + r 77^ - M r 77^ 



(29) 



The different coefficients of the expansion are explained 
in the following. The linear term is due to the elastic 
energy of grain boundaries. The zero temperature con- 
tribution 7o is the elastic energy per unit surface of a flat 
or smooth grain boundary that, in the limit of low angle 
grain boundaries, is given hy^ 



70 



C66^. exD 
2ttD ^ 27r6 ' 



(30) 



V. GRAIN BOUNDARY INDUCED MELTING 

The stability of crystalline ordering in a vortex lattice 
beyond the well known Bragg Glass regime is still a mat- 
ter of investigation. Experimental results suggest that an 
increase in temperature above a certain critical value Tm 



where the x > factor takes into account core interac- 
tion effects. The 7^ term, on the other hand, accounts 
for thermal fluctuations. Indicating the elastic Hamil- 
tonian in Eq. (221) as H = / ^|^(/)(k)X(k)X(-k) with 
(/)(k) = e|k|, and e = -Kh'^KjlD^ ^ the partition function 
of a thermally perturbed grain boundary over a surface 
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S is thus 



and the corresponding free energy per unit surface 



7t 



(31) 



(32) 



with P — {KbT)~^. The above term can be determined 
exphcitly calculating the logarithm of the partition func- 
tion as 



InZ 
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1 D 
— - arctan — 
13^ a. 



1 0,z 

— arctan — 
at D 



(33) 



where we have introduced a short wavelength cutoff 
27r/az to delimit the integration domain along the z axis. 

The r coefficient of the term is proportional to the 
energy of a junction between two grain boundaries and 
details of its computation will be given in Appendix. 

The term captures the case of the intermission of a 
third grain boundary in a junction, screening the effect 
introduced by the contribution. When this is the case, 
one loses an energy equal to Tn^ times the probability of 
such an event. In the low density limit, this probability 
is Mn, where M = 27r/D is roughly the interaction range 
of a grain boundarji2&. 

It is convenient to define Q = \nZ/S, so that the free 
energy functional in Eq. 129() can be rewritten as 



defining a melting temperature as T,„ = 



70 
KbQ' 



(34) 



As 



shown in Fig. |H1 for values of T close to T„i, f{n) 
shows a global minimum corresponding to a GB den- 

1 1 + y/l + iKBQiTn, - T)M/T 



sity 



where R is the 



R 3M 
average grain size. As discussed above, this suggests the 

possibility of a polycrystalline arrangement before the 
amorphous phase takes over. As soon as T reaches its 
melting value T^, the global minimum density becomes 
of the order of D~^, grains cannot be defined, and the 
system loses polycrystalline ordering in favor of a liquid- 
amorphous phase characterized by a typical dislocation 
spacing of order a. 

The considerations above allow to draw a phase dia- 
gram for the vortex array at low applied magnetic fields 
(i.e. when effects of disorder can be neglected). The re- 
sulting plot is shown in Figure The melting line is ob- 
tained plotting the above temperature T„i as a function 
of the magnetic induction. Here we use the expression 
for the local value of cge reported in Ref. I53 The curve 
shows reentrant behavior expected for low fields, due to 
the exponential decay of the elastic shear modulus in the 
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FIG. 6: Free energy density as a function of grain boundary 
density close to thermal melting point. 



B/Bc2 ^ limit. The line delimiting lattice and poly- 
crystal phases, instead, is obtained imposing that the free 
energy minimum shown in Fig.|Blis a global minimum. In 
the presence of disorder, we obviously expect modifica- 
tions of this schematic phase diagram. Nevertheless, for 
weak enough disorder the main features should remain 
valid. 
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FIG. 7: Phase diagram of the vortex ensemble for low values 
of the reduced field B/Bc^- The temperature is rescaled by 
the quantity a = Kb/Uso), where eo = ($o)^/(47rA)^ is an 
energy per unit length along the magnetic field direction, i.e. 
the typical energy for vortex interactions. The melting line is 
anticipated by the insurgence of a polycrystalline ordering. 
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VI. SUMMARY AND OUTLOOK 

In this paper we have investigated the properties of 
grain boundaries in a vortex polycrystaUine phase. A 
vortex polycrystal is experimentally observed in field 
cooling experiments when grain growth is arrested by 
disorderii^iiSiiLi2ii2i2S, and could also arise close to the 
melting line as an intermediate stage between the vortex 
lattice and the vortex liquid or glass-^'^*'. In both cases, 
the dynamics of the system can be studied analyzing 
grain boundaries, which play a similar role to that of do- 
main walls in ferromagnetic systems. Grain boundaries 
can be seen as elastic manifolds whose non-local surface 
tension can be obtained from the elastic description of 
the vortex lattice. Deformations are due to thermal fluc- 
tuations or to random stresses induced by vortex lattice 
deformationsSS. Once the main ingredients (i.e. elasticity 
and disorder) have been properly described, grain bound- 
aries can be studied with standard scaling methods, used 
in the past for various systems from flux lines to ferro- 
magnetic domain walls. In particular, we have studied 
disorder induced roughening, depinning under an applied 
stress and creep. These results are important to quantify 
arrested grain growth due to disorder and can be used to 
estimate the grain size in field cooling experiments^. 

An important question concerns the relevance of a 
polycrystaUine vortex structure for the transport proper- 
ties of a superconductor—. We have shown by numerical 
simulations that the critical current of a vortex polycrys- 
tal is systematically higher than the one observed in the 
corresponding single crystal case. This result reflects the 
fact that a polycrystal is pinned more effectively than a 
single crystal. 

Finally, we have extended the theory of grain bound- 
ary induced meltingSS to vortex lattices. We have written 
a free energy as a function of the grain boundary den- 
sity considering the contributions due to thermal fluctu- 
ations, elastic deformations, and junction formation. We 
find that the ordered crystal melts into a liquid passing 
through an intermediate polycrystaUine phase in agree- 
ment with recent numerical results'^^'^^. We have drawn 
a schematic phase diagram as a function of temperature 
and magnetic field which, however, can only be consid- 
ered as a first rough approximation. We have not taken 
into account the effect of disorder, which is believed to 
be responsible for a field induced transition to an amor- 
phous vortex glass. In addition, we have neglected the 
effect of isolated dislocations and their interactions with 
grain boundaries. Therefore at this stage the present the- 
ory should be seen mainly as a framework for a general 
physical mechanism^, supported by simulations'^^'^^ and 
by some experiments'*®, for vortex lattice melting. 
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Appendix 



ESTIMATE OF THE ENERGY ASSOCIATED TO 
JUNCTION FORMATION 

The presence of a term in the free energy func- 
tional [5^ was first suggested by Chui^® , in order to take 
into account grain boundary crossing in the framework 
of a crystal melting theory. Such a crossing energy con- 
sisted of a thermal contribution due to coupling between 
fluctuations of dislocations of crossing grain boundaries. 
Nonetheless, Bitter decoration experiments show that in 
vortex polycrystals, grain boundaries primarily rearrange 
forming junctions, instead of simply crossing. The forma- 
tion of such junctions determines variations in the overall 
free energy of the system due to two different contribu- 
tions, a zero temperature junction elastic energy and a 
thermal part related to fluctuations. In the following, we 
will address to these contributions respectively as Fq and 
Ft, being F = Fq -f Ft. 



1. Zero temperature energy 

We assume that because of the short range nature of a 
grain boundary stress field, grain boundary interactions 
are screened for long distances and we show that forming 
a junction leads to a zero-temperature elastic energy gain 
Fo^O. 

The idea is to focus on what happens when two grain 
boundaries come so close that they can form a junction. 
Let us consider a first grain boundary, e.g. directed along 
j with Burgers vectors b„ such that, b„ • i = 6„, and a 
single dislocation, belonging to the other grain boundary, 
whose Burgers vector is b' • i = —b'costp, being ip the 
junction angle. 

Since grain boundary interactions are short-ranged, we 
expect misorientations effects to make no difference in 
the energy computation until dislocations come close to 
a distance that we will call s. If, on a distance s, the 
interaction energy for (/j = is lower than for ^ 0, 
there in no reason for the system to make a junction. 
Otherwise, if there is an energy gain, grain boundaries 
are likely to join. 
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Considering the general expression for dislocation interactions 
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where r and r' are the positions of interacting dislocations, the energy (per unit length) of our system (GB and rotated 
dislocation) is 
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where s is the distance between the rotated dislocation and the grain boundary. Moreover, after summing the series. 
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Assuming that we have M dislocations within the range boundaries in a junction can also be estimated follow- 
ing Rcf. 38, After performing the thermal average of the 
interaction potential (|35|l . calculated on the cylinder of 



of s, the energy gain due to a junction will be 

M 



(38) 



radius |r| < (2/-\/3)-D/7r and taking the short range log- 
arithmic part of V 



Since the stress field generated by a grain boundary is 
exponentially suppressed beyond a distance of the same 
order of the dislocation spacing, we can give a rough es- 
timate of the sum taking M = I and si — D, i.e. 



Kb^ 
"27 
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(39) 



Thermal fluctuations 



(40) 



being Af a normalization constant. Evaluating the inte- 
gral for T ~ Tm leads to 
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where (p is the average junction angle. In the estimate 
The Ft contribution, due to the coupling between of the lattice-polycrystal crossover we have assumed ip ~ 
fluctuations of dislocations belonging to different grain 7r/3, as it is often observed in decoration experiments. 
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